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Abstract 

We present the numerical investigation of the fermionic two-body decays of squarks 
in the Minimal Super symmetric Standard Model with complex parameters. In the 
analysis we particularly take into account the cosmological bounds imposed by WMAP 
data. We find that the phase dependences of the decay widths of the third family 
squarks, as well as those of the first and second families, are quite significant which 
can provide viable probes of additional CP sources. We plot the CP phase depen- 
dences for each fermionic two-body decay channel of squarks q, L (i=l,2 , q=u,d;c,s;t,b) 
and speculate about the branching ratios and total (two-body) decay widths. 

PACS numbers: 14.80.Ly, 12.60. Jv 

1 Introduction 

The experimental HEP frontier is soon reaching TeV energies and most of physicists ex- 
pect that just there theoretically proposed Higgs bosons and super partners are waiting 
to be discovered. There are many reasons to be so optimistic. First of all, in spite of its 
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remarkable successes, the Standard Model has to be extended into a more complete the- 
ory which should solve the hierarchy problem and stabilize the Higgs boson mass against 
radiative corrections. The most attractive extension to realize these objectives is supersym- 
metry (SUSY) [1]. Its minimal version (MSSM) requires a non-standard Higgs sector [2] 
which introduces additional sources of CP- violation [3,4] beyond the Sckm phase [5]. The 
plethora of CP-phases also influences the decays and mixings of B mesons (as well as D 
and K mesons). The present experiments at BABAR, Tevatron and KEK and the one to 
start at the LHC will be able to measure various decay channels to determine if there are 
super symmetric sources of CP violation. In particular, CP-asymmetry and decay rate of 
B — > X s ^ form a good testing ground for low-energy supersymmetry with CP violation [6] . 
The above-mentioned additional CP-phases explain the cosmological baryon asymmetry of 
the universe and the lightest SUSY particle could be an excellent candidate for cold dark 
matter in the universe [7, 8] . 

In the case of exact supersymmetry, all scalar particles would have to have same masses 
with their associated SM partners. Since none of the superpartners has been discovered, 
supersymmetry must be broken. But in order to preserve the hierarchy problem solved the 
supersymmetry must be broken softly. This leads to a reasonable mass splittings between 
known particles and their superpartners, i.e. to the superpartners masses around 1 TeV. 

The precision experiments by Wilkinson Microwave Anisotropy Probe (WMAP) [9] 
have put the following constraint on the relic density of cold dark matter 3 

0.0945 < VLcDhih 2 < 0.1287 (1.1) 

Recently, in the light of this cosmological constraint an extensive analysis of the neu- 
tralino relic density in the presence of SUSY-CP phases has been given by Belanger et 
al. [10]. 

In this study we present the numerical investigation of the fermionic two-body decays 
of squarks in MSSM with complex SUSY parameters. Actually, we had performed a 



3 In our calculation, we have used WMAP-allowed bands in the plane Mi — ip which are based on 1st year 
data. Now the WMAP 3rd year data is also available [11], but the new WMAP + SDSS combined value for 
relic density of dark matter does not change the numerical results in Ref. [10], namely the WMAP-allowcd 
bands. See "Note added" section of Ref. [10]. 
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short study in this direction for the third family squarks [12] incorporating all the existing 
bounds on the SUSY parameter space by utilizing the study by Belanger et al. [10] . This 
investigation showed us that the effects of Mi and its phase (pu(i) on the decay widths of 
ti j2 and 6i j2 are quite significant. Now we extend it to all three families. Although the 
SUSY parameters fx, M 1: M 2 and Af are in general complex, we assume that fx, M 2 and 
Af are real, but Mi and its phase ipu(i) take values on the WMAP-allowed bands. These 
bands also satisfy the EDM bounds [13]. The experimental upper limits on the EDMs of 
electron, neutron and the 299 Hg and the 205 T/ atoms may impose constraints on the size 
of the SUSY CP-phases [14,15]. However, these constraints are highly model dependent. 
This means that it is possible to suppress the EDMs without requiring the various SUSY 
CP-phases be small. For example, in the MSSM assuming strong cancelations between 
different contributions [16], the phase of \x is restricted to |</? M | < n/10, but there is no such 
restriction on the phases of Mi and Af. In addition, we evaluate the parameter M 2 via 
the relation M 2 = (3/5)\M 1 \(t&n6 w )~ 2 which can be derived by assuming gaugino mass 
unification purely in the electroweak sector of MSSM. It is very important to insert the 
WMAP-allowed band in the plane Mi — (p into the numerical calculations instead of taking 
one fixed Mi value for all ^-phases, because, for example, on the allowed band for [i = 200 
GeV, Mi starts from 140 GeV for <p = and increasing monotonously it becomes 165 GeV 
for (p — 7r. In Ref. [10] two WMAP-allowed band plots are given, one for \i = 200 GeV and 
the other for /i = 350 GeV. For both plots the other parameters are fixed to be tan (3 = 10, 
m H + = 1 TeV, Af = 1.2 TeV, <ffj,=<fA f =0- We here simply choose generic masses for 
squarks as mg 2 =1000 GeV and 77^=750 GeV (q=u,d,c,s,t,b), although we had chosen 
more reasonable mass values in our previous note [12]. Taking different sets of values for 
squark masses around 1 TeV does not affect much neither the phase dependences nor the 
decay widths. As a second possibility, we take masses of 10 TeV for the first and second 
family squarks. 
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2 Squark Masses, Mixing and Two-body Decay Widths 
2.1 Masses and mixing in squark sector 

The superpartners of the SM fermions with left and right helicity are the left and right 
sfermions. In the case of top squark (stop) and bottom squark (sbottom) the left and right 
states are in general mixed. Therefore, the sfermion mass terms of the Lagrangian are 
described in the basis {qiAo) as [17, 18] 



M 2 RL M 2 RR 



with 



M\ L = M| + (7f L - e q sin 2 6 W ) cos(2[3)m 2 z + m 2 (2.2) 
M 2 RR = Mj, + e q sin 2 9 W cos{2(3)m 2 z + m 2 q (2.3) 
M 2 RL = (M 2 LR y = m q (A q -^ (tan f3)- 21 ^) (2.4) 

where m q , e q , l\ L and 6y/ are the mass, electric charge, weak isospin of the quarks and 
the weak mixing angle, respectively, tan/5 = v-ijv\ with V{ being the vacuum expectation 
values of the Higgs fields H®, i = 1,2. The soft-breaking parameters Mq, Mq, = M^(M^) 
for up components (down components), A q involved in Eqs. (2.2-2.4) can be evaluated for 
our numerical calculations using the following relations 



M l = + 



( sm 2 1 

v 2 3 

1 



2 



/(m| -m|) 2 - 


Am 2 q \A q 


- fi* cot/9| 2 ) 


cos(2/3)m 2 — m 2 






km\-m\) 2 - 


4m 2 q \A q 


-fjL* COt /?| 2 ) 


X)m\ - m\ 







(2.5) 



(2.6) 



M l = \ ( m l + m i =F \l ( m l - m l¥ - 4m t\ A i - V* cot/5| 2 ) 

+\ sin 2 6 W cos(2/?)m 2 - m\ (2.7) 

The squark mass eigenstates q\ and q\ can be obtained from the weak states qi and q R 
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via the g-mixing matrix 



e l(p v cos 6q sin 6q 



TV = . _■ . (2-8) 



— sin 9q e %Lp * cos 9 



where 



and 



<PS = arg[M 2 L ] = arg L4 ? - (2.9) 



cos 6q = sin 6q = LL qi (2.10) 

^I^P + K-^l) 2 ' yJ\Ml R \* + (rnl-Ml L )* 

One can easily get the following squark mass eigenvalues by diagonalizing the mass matrix 
in Eq. (2.1): 

< 2 = I ( M Il + m rr =F sJ{Ml L -Ml R f + A\Ml R \^ , m qi < m~ q2 (2.11) 

Note that the left-right mixing is not significant for the first and second generations; 
i.e. , qi ~ q R and q 2 ~ (q=u,d,c,s). 

We might add a comment about the possibility of a flavor mixing, for example, between 
the second and third squark families. In this case, the sfermion mass matrix in Eq. (2.1) 
becomes a 4x4 matrix in the basis (cl, c r , f^, i R ). Then one obtains squark mass eigenstates 
(ci, c 2 , ti, t 2 ) from these weak states, and analyzes their decays by utilizing procedures 
similar to the ones indicated in the text. The problem with flavor violation effects is that 
their inclusion necessarily correlates B, D and K physics with direct sparticle searches at 
colliders. Moreover, it has been shown that, with sizeable supersymmetric flavor violation, 
even the Higgs phenomenology at the LHC correlates with that of the rare processes [19]. 
In this work we have neglected such effects; however, we emphasize that inclusion of such 
effects can give important information on mechanism that breaks supersymmetry via decay 
products of squarks. 

2.2 Fermionic decay widths of squarks 

The quark- squark- chargino and quark-squark-neutralino Lagrangians have been first given 
in Ref. [1]. Here we use them in notations of Ref. 



= 9u{&Pr + k%P L )xUi + gd(l%P R + 4.P L )x+ c ^ + h.c. (2.12) 



and 



= 9*Wr + b\ k P L )xl<li + h.c (2.13) 

where u (u) stands for up-type (s)quark and d (d) stands for down-type s(quark). We 
also borrow the formulas for the partial decay widths of (ji ( qi=ti,Ui,bi and di) into quark- 
chargino (or neutralino) from Ref. [20]: 



g 2 \ l / 2 {ml,m 2 ,,m 2 ± ) 

r «' ~' + ® = Ifrrmj X 

I4I 2 + l4l 2 )K - - - 4ite(*£4)m,m^2.14) 



and 



g 2 X 1/2 {m 2 - m 2 m 2 ) 

r «.^ + *°> = iw x 

~(l4| 2 + l4| 2 )K ~m 2 q - m|o) - 4i^4*4K™x°] (2.15) 

with A(x, y, = x 2 + y 2 + z 2 — 2{xy + xz + yz). 
The explicit forms of £f k , k q ik are 

4 = -<^i + W 2 V fe2 , kl = K«lY D U* k2 , (2.16) 
for up squarks u, c, t and 

4 = -n*^i+^4*£4 2 , 4 = n*W 2 , (2.i7) 

for down squarks <i, s, 6 (U=w, c, t ; D=<i, s, b). And a? fc , are given as 

4 = E(^L)^L, 4 = X>l)*£L (2.i8) 

n=l n=l 

where 

^ = (©- < 219 » 

Here, the components for up squarks u, c, t are 

fl k = -^(N k2 + ^t a n9 w N kl ) 
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(2.20) 



and for down squarks d, s, b 




1 



(N k2 



tan 9 w N kl ) 



V2 




— — tan 6 W N^ 
(h g Rk Y = ~Y D N* k: 



hi 



Lk 



(2.21) 



We have to point out that although at the loop level the SUSY-QCD corrections could 
be important, our analysis here are merely at tree level, as can be seen from Eqs. (2.13) 
and (2.14). In this work we content with tree- level amplitudes as we aim at determining 
the phase-sensitivities of the decay rates, mainly. 

3 Two-body decays of first and second generations 
squarks 

Now we present the dependences of the u\, u 2 , di and d 2 two-body decay widths on the 
phase ipu(i) for li = 200 GeV and for /i = 350 GeV. Here we choose the values for the 
masses (m„ 2 = mj 2 , m„ 1 = ra-±, «^±, m x°) — (1000 GeV, 750 GeV, 180 GeV, 336 
GeV, 150 GeV) for li = 200 GeV and (m^ 2 = m^, = m ^-, m x*' m x?) = (1000 
GeV, 750 GeV, 340 GeV, 680 GeV, 290 GeV) for li = 350 GeV. Note that although the 
neutralino and chargino masses vary with (fiu(i), these variations are not large. Therefore, 
as a final state particle (i.e., on mass-shell), we have chosen fixed (average) mass values 
for charginos and neutralinos. Fig. 1(a) and Fig. 1(b) show the partial decay widths of the 
channels u 2 — * dxf, u 2 — > dx 2 , u 2 — > uxl, Ui — * ux± and d 2 — > uxi, d 2 — > ux 2 , d 2 — * dxl, 
d\ — > dxi, respectively for \x = 200 GeV. [ We plot the same processes for \x = 350 GeV in 
Fig. 2(a) and Fig.2(b) ]. The decay channels u\ — > dxi an d ^i ~~ * dx 2 are absent because 
the mixing elements (TZij) and the Yukawa couplings Y u , Yd in £fj and kfj parameters are 
zero, whereas the channel U\ uxl * s present, because the parameter 6 n is not zero. 
Because of a similar reason the phase dependence of Hi — > uxl channel is very little. 
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The branching ratios for u 2 are roughly B(u 2 — > c^X^) : S(w 2 — > dxt) '■ B{u 2 — > MX?) 
^5:2: 0.05. for the case \x = 200 GeV. From Fig. 1(b) one can see that B(d 2 — > u\ 2 ) 
: B(d 2 -> wxr) : #(4 -> ~ 6 : 1 : 0.2 for the case fj, = 200 GeV. The width 

r(u 2 — > dx 2 ) increases as the phase increases from to n, but T(u 2 — > dxt) decreases as 
<fu(i) increases, both showing significant dependence on the phase. The phase dependence 
is more significant for the decay channel u 2 — > uxl (its decay width increases from 0.0006 
GeV to 0.15 GeV as <fu(i) increases from to it) but its branching ratio is too small to 
observe it. 

The channels d\ — > uxi and d\ — > ux 2 are absent because the related mixing elements 
and Yukawa couplings are zero ( therefore £f ± , kf ± , £f 2 , kf 2 become zero); whereas the 
channel d x — > dx\ is present (since bf^O), but it is little with respect to d 2 decay channels. 
Its phase dependence is also not so strong ( from (fiu(i)= to tt, its width increases only 
from 0.264 GeV to 0.271 GeV ). Therefore, the d\ decay can be ignored in comparison with 
the d 2 decay. 

The decay pattern of the second family squarks coincides with that of the first sfamily, 
because we can set m s m0 and m c m0 like m u m0 and m^O in comparison with the masses 
mg 2 =1000 GeV, 77^=750 GeV ( The use of m s ^0.1 GeV and m c ~1.5 GeV does change 
nothing). 

In the analysis of cosmological constraints by Belanger et al. they take masses of 10 TeV 
for the first and second generations of squarks. As an academic exercise, we have repeated 
the above calculations only changing the first and second squarks masses to m„ 2 = mj 2 = 10 
TeV and = TeV. We obtained again similar phase dependences, but larger 

(~ one order) width values because of huge phase spaces. 

4 Two-body decays of third generation squarks 

In this section we plot the dependences of the t±, t 2 , b\ and b 2 partial decay widths on tputi) 
for n = 200 GeV and for \x = 350 GeV. Here we simply choose mass values for the third 
family squarks as mi 2 =mi 2 =1000 GeV and =777^=750 GeV in order to compare the 
results with those of the first two families. We had taken more reasonable mass spectrum 
for squarks in our previous study [12]. The other masses are (^±, m x?) = 
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GeV, 336 GeV, 150 GeV) for /i = 200 GeV and (m~±, m~±, m^o) = (340 GeV, 680 GeV, 
290 GeV) for // = 350 GeV. 

For both sets of values by calculating the Mq and M v values corresponding to and 
mi 2 , we plot the decay widths for Mq > Mfj and Mq < Mfj, separately. Fig.3(a) and 
Fig. 3(b) show the partial decay widths of the channels ti — > bxt, h — > &x 2 ii — > 
?2 - ► &xf , ^2 - ► and ^2 - ¥ txi as a function of <fu{i) f° r — 200 GeV assuming 
Mq > My and Mq < M V) respectively. In these figures significant dependences on tpu(i) 
phase are seen. T(ti — > txt) and T(ti — > 6xi") decay widths increase as <fu(i) increases from 
to 7r, but r(^ — > 6%^) width decreases as ty'rm) increases. On the other hand, T(t 2 — > ixi) 
decrease for both Mq > Mfj and Mq < Mfj cases as ipu(i) increases; T(t 2 — > bxt) decreases 
for Mq > Mfy but increases for Mq < My and T(t 2 — ► &X2") increases for Mq > but 
decreases for Mq < M^-. 

The branching ratios for t 2 are roughly B(t 2 — >■ 6X1") : B(t 2 ^ txi) '■ B(t 2 — > 6%^) ~ 
12 : 2.5 : 1. This simply reflects both the large phase space and large Yukawa coupling for 
the decay t 2 — > bxt- 

In Figs. 4(a) and 4(b) you see the same partial decay widths for \x = 350 GeV. They, 
too, show the significant dependences on CP-violation phase. For both Mq > M v and 
Mq < M v cases; F(ti — > txt) and F(ti — > bxt) decay widths increase as (pu(i) increases 
from to 7r, but T(tx — > 6%^) width decreases as (fiu(i) increases. Besides that T(t 2 — > t%?) 
and F(t 2 — > 6%^ ) widths decrease but T(t 2 — ► &X2 ) width increases as (pu(i) increases from 
to 7r for both Mq > M v and Mq < M v cases. The branching ratios for t 2 are roughly 
B(t 2 - bxt) ■ B(t 2 - txi) ■ B(t 2 - &x+) » 8 : 2 : 1. 

For /i = 350 GeV the WMAP-allowed band [10] takes place in larger M\ values (~ 
305 — 325 GeV) leading to larger chargino and neutralino masses. This naturally leads 
very small decay width for t\ — > bxt ■ The decay width of the process t 2 — > bxt is the 
largest one among the t 2 channels and the branching ratios are B(t 2 — > bxt) '■ B(t 2 — > tx?) 
: _B(t 2 — > 6% J) « 8 : 3 : 1. The decay i 2 — > shows strong phase dependence. 

We give sbottom decay widths as a function of (fiu(i) in Figs. 5 (a)- (b) (for \x = 200 GeV 
) and in Figs. 6(a)-(b)(for /j, = 350 GeV ). T(b 2 — > 6x?) is smaller than T(6 2 — * tx~) i n 
spite of large phase space, because in b 2 — > 6x5 decay only Y& coupling enters which is very 
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small in comparison with Y t . The dependences of the phase ipu(i) in sbottom decays are 
similar to those in stop decays. 

For the case Mq > Mp, b 2 decay is 7-10 times larger than that of 61, whereas for 
Mq < Mi) the reverse is true. The branching ratios for b 2 decays are B{b 2 — > tx±) '■ 
B(b 2 -> : B(b 2 -> ~ 14 : 10 : 0.5 for fi = 200 GeV and Mq > M 3 , and similarly 
14 : 3 : 0.5 for /i = 350 GeV and Mq > M^. While the process b 2 — > bx\ is suppressed 
more than one order its dependence on (fiu(i) is prominent such that the value of decay 
width at (fu(i) = is 2 times larger than that at <fu{i) = (^-dependences of the processes 
i>2 — * tXi and b 2 — > t\ 2 can be seen easily in Fig. 5(a). 

The branching ratios for b\ decays are B(b\ — > tx±) '■ B(b\ — > tx 2 ) '■ B(bi — > bx^) ~ 9 
: 3 : 0.5 for fx = 200 GeV and Mq < M b and -> t^r) : ~ 7 : 0.3 for 

^ = 350 GeV and M Q < M D . 

5 Discussion 

In this paper, we present the numerical investigation of the fermionic two-body decays of 
squarks in three families in the Minimal Supersymmetric Standard Model with complex 
parameters taking into account the cosmological bounds imposed by WMAP data. Nu- 
merical calculations of decay widths at tree level show significant dependence on the CP 
phase tputi) for the third family squarks, as well as for the first and second families. 

We have assumed the same mass of 750 GeV, for all squarks of type 1 and the same 
mass of 1000 GeV for all squarks of type 2. The decay width values of the third family 
squarks and those of the first two families differ only two to three times in favor of the 
third family. 

In the case of \x = 200 GeV all the channels of u 2 decay (d 2 decay), i.e. , u 2 — > dxt, 
u 2 — > dxf and u 2 — > uxi (d 2 — > ux 2 , d 2 — > UX\ and d 2 — > dxl) are present and have very 
significant phase dependences, but the channel u 2 — > ux\ (d 2 — > dx\) has a very small 
branching ratio. At the u\ decay {d\ decay) only one channel, i.e., u\ — > uxi (di — * dxi) 
is open with a small branching ratio. The decay channels u\ — > dxt and u\ — > dxt 
(di — > uxi and d 1 — * ux 2 ) are absent because of the reason explained in section 3. 

Very roughly speaking, for /j, = 200 GeV and Mq > Mfj or Mp we find the following 
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total (two-body) widths: 



r t otai(w2orc 2 ) w T totBl (d 2 OTS2) ^7GeV (5.1) 

rtatai(wiorci) ~ lGeV (5.2) 

iWdiorsO ~ 0.5 GeV (5.3) 

Ttotaife) » 15 GeV (5.4) 

r t otai(6 2 ) « 24 GeV (5.5) 

r to tai(ti) « 10 GeV (5.6) 

r to tai(&i) « 3 GeV (5.7) 



From these results we see that although the decays of the third family squarks are 
more important, the decays of the first two families are not ignorable at all. For example, 
if we assume that the probability of producing every kind of squarks in a proton-proton 
collision (LHC) are more or less equal, then the ratio of total decay width of the third 
family squarks to that of the first two families would be approximately 2. 

In the case of the choice of squarks masses of 10 TeV for the first two families and 1 
TeV for the third family which we call it as "an academic exercise" (because there is no 
hope to produce the first two family squarks in LHC) this ratio even reverses in favor of 
the first two sfamilies. 
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Figure 1: (a)-(b) Partial decay widths T of the Ui^ and c? 12 decays for fi = 200 GeV , 
tan/5 = 10, A u =A d =1.2 TeV, ip^=ip Au =ip Ad =0, 111^=111^ = 750 GeV and mr a ,=m^= 1000 
GeV. 
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Figure 5: (a)-(b) Partial decay widths T of the &i j2 decays for fi = 200 GeV , tan/? 
A 6 =1.2 TeV, ip^=(p Ab =0, m 5l = 750 GeV and mgj = 1000 GeV. 
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